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The Stochastic Calculus of Looping Sequences (SCLS) is a recently proposed modelling language for 
the representation and simulation of biological systems behaviour. It has been designed with the aim 
of combining the simplicity of notation of rewrite systems with the advantage of compositionality. It 
also allows a rather simple and accurate description of biological membranes and their interactions 
with the environment. 

In this work we apply SCLS to model a newly discovered anmionium transporter. This trans- 
porter is believed to play a fundamental role for plant mineral acquisition, which takes place in the 
arbuscular mycorrhiza, the most wide-spread plant-fungus symbiosis on earth. Due to its potential 
application in agriculture this kind of symbiosis is one of the main focuses of the BioBITs project. 

In our experiments the passage of NH^ I NH^ from the fungus to the plant has been dissected 
in known and hypothetical mechanisms; with the model so far we have been able to simulate the 
behaviour of the system under different conditions. Our simulations confirmed some of the latest 
experimental results about the LjAMT2;2 transporter. The initial simulation results of the modelling 
of the symbiosis process are promising and indicate new directions for biological investigations. 

1 Introduction 

Given the central role of agriculture in worldwide economy, several ways to optimize the use of costly 
artificial fertilizers are now being actively pursued. One approach is to find methods to nurture plants 
in more "natural" manners, avoiding the complex chemical production processes used today. In the last 
decade the Arbuscular Mycorrhiza (AM), the most widespread symbiosis between plants and fungi, got 
into the focus of research because of its potential as a natural plant fertilizer. Briefly, fungi help plants to 
acquire nutrients as phosphorus (P) and nitrogen (N) from the soil whereas the plant supplies the fungus 
with energy in form of carbohydrates [32]. The exchange of these nutrients is supposed to occur mainly 
at the eponymous arbuscules, a specialized fungal structure formed inside the cells of the plant root. The 
arbuscules are characterized by a juxtaposition of a fungal and a plant cell membrane where a very active 
interchange of nutrients is facilitated by several membrane transporters. These transporters are surface 
proteins that facilitate membrane crossing of molecules which, because of their inherent chemical nature, 
are not freely diffusible. 

As long as almost each cell in the majority of multicellular organisms shares the same genome, 
modem theories point out that morphological and functional differences between them are mainly driven 
by different genes expression [2]. Thanks to the last experimental novelties [40, 30] a precise analysis 
of which genes are expressed in a single tissue is possible; therefore is possible to identify genes that 
are pivotal in specific compartments and then study their biological function. Following this route a new 
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membrane transporter has been discovered by expression analysis and further characterized [18]. This 
transporter is situated on the plant cell membrane which is directly opposed to the fungal membrane, 
located in the arbuscules. Various experimental evidence points out that this transporter binds to an NH^ 
moiety outside the plant cell, deprotonates it, and mediates inner transfer of NH^, which is then used 
as a nitrogen source, leaving an //+ ion outside. The AM symbiosis is far from being unraveled: the 
majority of fungal transporters and many of the chemical gradients and energetic drives of the symbiotic 
interchanges are unknown. Therefore, a valuable task would be to model in silico these conditions and 
run simulations against the experimental evidence available so far about this transporter. Conceivably, 
this approach will provide biologists with working hypotheses and conceptual frameworks for future 
biological validation. 

In computer science, several formalisms have been proposed for the description of the behaviour of 
biological systems. Automata-based models [3, 28] have the advantage of allowing the direct use of 
many verification tools such as model checkers. Rewrite systems [11, 35, 8, 1] usually allow describing 
biological systems with a notation that can be easily understood by biologists. Compositionality allows 
studying the behaviour of a system componentwise. Both automata-like models and rewrite systems 
present, in general, problems from the point of view of compositionality, which, instead, is in general 
ensured by process calculi, included those commonly used to describe biological systems [36, 34, 9]. 

The Stochastic Calculus of Looping Sequences (SCLS) [6] (see also [8, 7, 29]) is a recently proposed 
modelling language for the representation and simulation of biological systems behaviour. It has been 
designed with the aim of combining the simplicity of notation of rewrite systems with the advantage of a 
form of compositionality. It also allows a rather simple and accurate description of biological membranes 
and their interactions with the environment. 

In this work we apply SCLS to model the mentioned transporter. This transporter is differentially 
expressed in arbuscular cells and is believed to play a fundamental role in the nutrients uptake which 
takes place in the context of plants-fungi symbiosis at the root level. This symbiosis is one of the main 
focuses of the BioBITs project, due to its relevant role in agriculture. 

On these premises, the aim of this work is to model the interchange between the fungus-plant inter- 
face and the plant cells, using as a reference system the NH^/NH^ turnover. This could disclose which 
is the driving power of the net nitrogen flux inside plant cells (still unknown at the chemical level). 
Furthermore, this information may also be exploited to model other, and so far poorly characterized, 
transporters. 

Outline Section 2 introduces the syntax and the semantics of the SCLS and shows some modelling 
guidelines. Section 3 presents the SCLS representations of the ammonium transporter and our experi- 
mental results, discussed in Section 4. Section 5 concludes by outlining some further work. 

2 The Stochastic Calculus of Looping Sequences 

In this section we briefly recall the Stochastic Calculus of Looping Sequences (SCLS) [6]. A SCLS 
(biological) model consists of a term, representing the structure of the modelled system and a set of 
stochastic rewrite rules which determine its evolution. 

Terms Terms of SCLS are defined from a set of atomic symbols representing the basic elements of the 
system under description, which could correspond, according to the level of the description, to genes, 
proteins, molecules or even single chemical elements. Terms are built from the atomic elements via 
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the operators of sequencing (defining an ordered sequence of atomic elements), looping (defining mem- 
branes from looping sequences) and parallel composition (representing coexistence in the same biologi- 
cal ambient). 

Assuming a possibly infinite alphabet S' of symbols ranged over by a^b^c^... the syntax of SCLS 
terms T and sequences S is defined by the grammar: 

T ::= S \ {S)^\T \ T\T 
S ::= e \ a \ S-S 

where a is any element of co and £ is the empty sequence. We denote the infinite sets of terms and 
sequences with ^ and respectively. 

The SCLS terms are built from the elements of S' by means of four syntactic operators: sequencing 
_ • _, looping (_)^, parallel composition _ | _, and the containment operator _ J _. Sequencing is used to con- 
catenate atomic elements. The looping operator transforms a sequence in a closed loop when combined 
with another term via the _J _ operator. A term can then be either a sequence, or a looping sequence 
containing another term, or the parallel composition of two terms. By the definition of terms, we have 
that looping and containment are always associated, hence we can consider (_)^J _ as a single binary 
operator that applies to one sequence and one term. 

The biological interpretation of the operators is the following: the main entities which occur in 
cells are DNA and RNA strands, proteins, membranes, and other macro-molecules. DNA strands (and 
similarly RNA strands) are sequences of nucleic acids, but they can be seen also at a higher level of 
abstraction as sequences of genes. Proteins are sequences of amino acids which typically have a very 
complex three-dimensional structure. In a protein there are usually (relatively) few subsequences, called 
domains, which actually are able to interact with other entities by means of chemical reactions. SCLS 
sequences can model DNA/RNA strands and proteins by describing each gene or each domain with a 
symbol of the alphabet. Membranes are closed surfaces often interspersed with proteins, and may have 
a content. A closed surface can be modelled by a looping sequence. The elements (or the subsequences) 
of the looping sequence may represent the proteins on the membrane, and by the containment operator 
it is possible to specify what the membrane contains. Other macro-molecules can be modeled as single 
alphabet symbols, or as sequences of their components. Finally, juxtaposition of entities can be described 
by the parallel composition operator of their representations. More detailed description of the biological 
interpretation of the SCLS operators involved in the description of the system presented in this paper will 
be given at the end of this section. 

Brackets can be used to indicate the order of application of the operators. We assume that (-)^J - 
has the precedence over _ | _ and _ • _ over _ | _, therefore (5)^ J Ti \ T has to be read as ((5)^ J Ti)\T and 
a • Z? I c as {a-b)\c. An example of an SCLS term is (3 1 | (m • n)^ \{c'd\e) consisting of three entities 
a, b and (m • n)^ \ {c- d\e). It represents a membrane with two molecules m and n (for instance, two 
proteins) on its surface, and containing a sequence c • d and a molecule e. Molecules a and b are outside 
the membrane. See Figure 1 for some graphical representations. 



Structural congruence As usual in systems of this kind, SCLS terms representing the same entity 
may have syntactically different structures. Such terms are thus identified by means of a structural 
congruence relation. In particular the structural congruence relations =s and =t are defined as the least 
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Figure 1: (i) (a-b-c)^ \ e; (ii) (a-b-c)^ \ {d-ef \ e; (iii) {a-b-cf\ {{{d-ef \ £)\f-g). 

congruence relations on sequences and on terms, respectively, satisfying the following rules: 

SiiSrS^) =s {SvS2)-S3 S£=s£S=sS 
Si =s S2 implies Si =t S2 and J T =t (52)^ J T 
Ti I T2 =T T2 I Ti Ti I {T2 I T3) =T {Ti I T2) \T3 T\£=tT 
Ti =T T2 implies (5)^ J Ti =t (5)^ J T2 
(e)^J £=t£ (5r52f J T =t (52-5if J T 

Structural congruence states the associativity of • and | , the commutativity of the latter and the neutral 
role of £. Moreover, axiom (5r52)^J T =t (S2-Si)^ \ T says that looping sequences can rotate. In the 
following we will simply use = in place of =t, =8- 

Rewrite Rules, Variables and Patterns A rewrite rule is defined as a pair of terms (possibly con- 
taining variables), which represent the patterns defining the system transformations, together with a rate 
representing the speed of the modelled reaction. Rules are applicable to all subterms, identified by the 
notion of reduction context introduced below, which match the left-hand side of the rule via a proper 
instantiation of its variables. The system transformation is obtained by replacing the reduced subterm by 
the corresponding instance of the right-hand side of the rule. 

Variables in patterns can be of three kinds: two are associated with the two different syntactic cate- 
gories of terms and sequences, and one is associated with single alphabet elements. We assume a set of 
term variables TV ranged over by X , F, Z, . . . , a set of sequence variables SV ranged over by x, z , . . . , and 
a set of element variables ^ ranged over by x,3;,z, . . .. All these sets are pairwise disjoint and possibly 
infinite. We denote by y the set of all variables TV USVU and with p any variable in y. A pattern 
is a term which may include variables, i.e. patterns P and sequence patterns SP of SCLS are given by 
the following grammar: 

P ::= SP I {SPf\P I P\P I X 
SP £ \ a \ SP'SP I 3c I X 

where a is an element of S', and X, x and x are elements of TV, SV and ^ , respectively. The infinite set of 
patterns is denoted with The structural congruence relation can be be trivially extended to patterns. 
An instantiation is a partial function a '.Y ^ ^ which preserves the type of variables, thus for X G 
ry,x G SV andx G 3^ we must have o{X) G ^, g{x) G and g{x) G respectively. Given P G ^, the 
expression Po denotes the term obtained by replacing each occurrence of each variable p G 1^ appearing 
in P with the corresponding term 0"(p). 
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Let Z denote the set of all the possible instantiations and Var{P) the set of variables appearing in 

P G Then a rewrite rule is a triple {Pi,P2,k), denoted with Pi Ap2, where k G IR-^ is the kinetic 
constant of the modeled chemical reaction, Pi,P2 G Pi ^ 8 such that Var{P2) C Var{Pi). A rewrite 

k 

rule Pi \-^P2 then states that a subterm Pi O", obtained by instantiating variables in Pi by some instantiation 
function a, can be transformed into the subterm P2G with kinetics given by the constant k. 

Contexts The definition of reduction for SCLS systems is completed by resorting to the notion of re- 
duction context. To this aim, as usual, the syntax of terms is enriched with a new element □ representing 
a hole. Reduction context (ranged over by C) are defined by: 

C ::= □ I C\T \ T\C \ {Sf \C 

where P G ^ and 5 G The context □ is called the empty context. We denote with the infinite set 
of contexts. 

By definition, every context contains a single □. Let us assume C^C G ^. Then C[T] denotes the 
term obtained by replacing □ with T in C; similarly C[C'] denotes context composition, whose result is 
the context obtained by replacing □ with C' in C. The structural equivalence is extended to contexts in 
the natural way (i.e. by considering □ as a new and unique symbol of the alphabet S). 

Note that the general form of rewrite rules does not permit to have sequences as contexts. A rewrite 

k 

rule introducing a parallel composition on the right hand side (as a Z? | c) applied to an element of a 
sequence (e.g., m-a-m) would result into a syntactically incorrect term (in this case m- {b\c) • m); which 
moreover does not seem to have any biological meaning, parallel and sequence represent orthogonal 
notions. Therefore to modify a sequence, a pattern representing the whole sequence must appear in the 

^ k I ~ 

rule. For example, rule a'X\-^a\x can be applied to any sequence starting with element a, and, hence, the 
term a-b can be rewritten as a\b, and the term a-b-c can be rewritten sls a\b'C. The notion of reduction 
context will forbid that a substitution like this is applied to looping sequences. 

Stochastic Reduction Semantics The operational semantics of SCLS is defined by incorporating a 
collision-based stochastic framework along the line of the one presented by Gillespie in [13], which is, 
de facto, the standard way to model quantitative aspects of biological systems. Following the law of mass 
action, it is necessary to count the number of reactants that are present in a system in order to compute the 
exact rate of a reaction. The same approach has been applied, for instance, to the stochastic ;r-calculus 
[33, 34]. The idea of Gillespie's algorithm is that a rate constant is associated with each considered 
chemical reaction. Such a constant is obtained by multiplying the kinetic constant of the reaction by the 
number of possible combinations of reactants that may occur in the system. The resulting rate is then 
used as the parameter of an exponential distribution modelling the time spent between two occurrences 
of the considered chemical reaction. 

The use of exponential distributions to represent the (stochastic) time spent between two occurrences 
of chemical reactions allows describing the system as a Continuous Time Markov Chain (CTMC), and 
consequently allows verifying properties of the described system analytically and by means of stochastic 
model checkers. 

The number of reactants in a reaction represented by a reduction rule is evaluated considering the 
number of occurrences of subterms to which the rule can be applied and of the terms produced. For 

, k 

instance in evaluating the application rate of the rewrite rule P = a | c to the term T = a\a\b\bv^Q 
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must consider the number of the possible combinations of reactants of the form a\b in T. Since each 
occurrence of a can react with each occurrence of b, this number is 4. So the appUcation rate of Risk -4. 

The evaluation of the application rate of a reduction rule containing variables is more complicate 
since there can be many different ways in which variables can be instantiated to match the subterm to be 
reduced, and this must be considered to correctly evaluate the application rate. The technique to do this 
is described in [6]. With this technique, given two terms T, and a reduction rule R, we can compute 
the function ff{R^T^T'), defining the number of the possible applications of the rule R to the term T 
resulting in the term ^^ We refer to [6] for more details and explanation. 

Given a finite set M of stochastic rewrite rules, then, the reduction semantics of SCLS is the least 
labelled transition relation satisfying the following rule: 

R = P^^P^^^ T = C[Pig] r = C[P2g] 

The rate of the reduction is then obtained as the product of the rewrite rate constant and the number of 
occurrences of the rule within the starting term (thus counting the exact number of reactants to which 
the rule can be applied and which produce the same result). The rate associated with each transition in 
the stochastic reduction semantics is the parameter of an exponential distribution that characterizes the 
stochastic behaviour of the activity corresponding to the applied rewrite rule. The stochastic semantics is 
essentially a Continuous Time Markov Chain (CTMC). A standard simulation procedure that corresponds 
to Gillespie's simulation algorithm [13] can be followed. The most recent implementation of SCLS, 
based on Gillespie's algorithm, is described in [38]. 

Modelling Guidelines SCLS can be used to model biomolecular systems analogously to what is done, 
e.g, by Regev and Shapiro in [37] for the ;r-calculus. An abstraction is a mapping from a real-world 
domain to a mathematical domain, which may allow highlighting some essential properties of a system 
while ignoring other, complicating, ones. In [37], Regev and Shapiro show how to abstract biomolec- 
ular systems as concurrent computations by identifying the biomolecular entities and events of interest 
and by associating them with concepts of concurrent computations such as concurrent processes and 
conmiunications. 

The use of rewrite systems, such as SCLS, to describe biological systems is founded on a different 
abstraction. Usually, entities (and their structures) are abstracted by terms of the rewrite system, and 
events by rewrite rules. We have already recalled the biological interpretation of SCLS operators in the 
previous section. 

In order to describe cells, it is quite natural to consider molecular populations and membranes. 
Molecular populations are groups of molecules that are in the same compartment of the cell. As we 
have said before, molecules can be of many types: they could be classified as DNA and RNA strands, 
proteins, and other molecules. Membranes are considered as elementary objects, in the sense that we 
do not describe them at the level of the lipids they are made of. The only interesting properties of a 
membrane are that it may have a content (hence, create a compartment) and that it may have molecules 
on its surface. 

We give now some examples of biomolecular events of interest and their description in SCLS. The 
simplest kind of event is the change of state of an elementary object. Then, we consider interactions 
between molecules: in particular complexation, decomplexation and catalysis. These interactions may 
involve single elements of non-elementary molecules (DNA and RNA strands, and proteins). Moreover, 



M. Coppo, F. Damiani, E. Grassi, M. Guether and A. Troina 



83 



Biomolecular Event 


Examples of SCLS Rewrite Rule 


State change 


a b 

x-a-y X'b-y 


Complexation 


a\b ^ c 

x-a-ylb i-> X'C-y 


Decomplexation 


c a\b 

x-c-y x-a-y\b 


Catalysis 


c\Pi c\P2 (where Pi P2 is the catalyzed event) 


Membrane crossing 


a\ (3f)'^JX ^ (3f)'^J(a|X) 
{5^^\{a\X) ^ a\ (3f)^JX 
x-a-y\ {I)^\X ^ (^^\{x-a-y\X) 
{2)^\(I-a-y\X) ^ x-a-y\ (?)^JX 


Catalyzed 
membrane crossing 


a\ {b'^'^lX ^ (Z?-3f)'^J(a|X) 
{b'I)^\{a\X) ^ a\ ib'^^\X 
x-a-y\ {b-z)^\X ^ {b -z)^ \ (x- a-y\X) 
{b-I)'^\ {x-a-y\X) ^ x-a-y\ ib-^^\X 


Membrane joining 


i^'^\ia\X) ^ (a-^^JX 
i5^^\{y-a-z\X) ^ {y-a-z-^^\X 


Catalyzed 





Table 1: Guidelines for the abstraction of biomolecular events into SCLS. 



there can be interactions between membranes and molecules: in particular a molecule may cross or join 
a membrane. 

Table 1 lists the guidelines (taken from [6]) for the abstraction into SCLS rules of the biomolecular 
events we will use in our application.^ Entities are associated with SCLS terms: elementary objects are 
modelled as alphabet symbols, non-elementary objects as SCLS sequences and membranes as looping 
sequences. Biomolecular events are associated with SCLS rewrite rules. 

3 Modelling the Ammonium Transporter with SCLS 

The scheme in Figure 2 (taken from [18]) illustrates nitrogen, phosphorus and carbohydrate exchanges at 
the mycorrhizal interface according to previous works and the results of [18]. In this paper we focus our 
investigation on the sectors labelled with (c), (al) and (a2). Namely, we will present SCLS models for the 
equilibrium between NH^ and NH3 and the uptake by the LjAMT2;2 transporter (c), and the exchange 
of NH^ from the fungus to the interspatial level (al-2). The choice of SCLS is motivated by the fact that 
membranes, membrane elements (like LjAMT2;2) and the involved reactions can be represented in it in 
a quite natural way. 

The simulations illustrated in this section are done with the SCLSm prototype simulator [38]. In 
the following we will use a more compact notation for the parallel composition operator, namely, we 
will write axn to denote the parallel composition of n atomic elements a. Moreover, to simplify the 
counting mechanism, we might use different names for the same molecule when it belongs to different 



Kinetics are omitted from the table for simplicity. 
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Non-Arbuscule Primary and Arbuscule containing 
containing root secondary root cortex cell 

cortex cell cell wall 



Interfacial apoplast/ 
Periarbuscular Space 



Arbuscule branch 




(g) 



Sucrose 

Pi 



- Sucrose Hexose 

Pi 



Hexose 



, (d) 



(e) — H+p- — (f) 



Phosphate 




Phosphate 


transporter 






transporter 


(plant) 






(fungus) 



Pi 

(f) ■ 



(b) 



NH, 



NH, 



H+-ATPase 
(plant) 



H+-ATPase 
(fungus) 




I NH 



3 : . 



(b) 



NH, 



(a1) 



(c) 



AA * 



NH, 



. NH, 



LjAMT2;2 



NH4" 



-(^^^^^^ (a2) 



1 



Plant membranes 



Plant membrane/ periarbuscular Fungal membrane 
membrane (PAM) 




Figure 2: (al-2) NH^/NH^ is released in the arbuscules from arginine which is transported from the extra- to the in- 
traradical fungal structures [16]. NH^/NH^ is released by so far unknown mechanisms (transporter, diffusion (al) or 
vesicle-mediated (a2)) into the periarbuscular space (PAS) where, due to the acidic environment, its ratio shifts towards 
NH^ (> 99.99%). (b) The acidity of the interfacial apoplast is established by plant and fungal //+-ATPases [21, 5] thus 
providing the energy for -dependent transport processes, (c) The NH^ ion is deprotonated prior to its transport across the 
plant membrane via the LjAMT2;2 protein and released in its uncharged NH^ form into the plant cytoplasm. The NH^/NH^ 
acquired by the plant is either transported into adjacent cells or immediately incorporated into amino acids (AA). (d) Phos- 
phate is released by so far unknown transporters into the interfacial apoplast. (e) The uptake of phosphate on the plant side 
then is mediated by mycorrhiza- specific Pi-transporters [25, 18]. (f) AM fungi might control the net Pi-release by their own 
Pi-transporters which may reacquire phosphate from the periarbuscular space [5]. (g) Plant derived carbon is released into 
the PAS probably as sucrose and then cleaved into hexoses by sucrose synthases [22] or invertases [39]. AM fungi then ac- 
quire hexoses [41, 42] and transport them over their membrane by so far unknown hexose transporters. It is likely that these 
transporters are proton co-transporter as the GpMSTl described for the glomeromycotan fungus Geosiphon pyriformis [31]. 
Exchange of nutrients between arbusculated cells and non-colonized cortical cells can occur by apoplastic (h) or symplastic 
(i) ways. 
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compartments. For instance, occurrences of the molecule NH3 inside the plant cell are called NH3 Jnside. 
3.1 NH^/NH^ Equilibrium 

We decided to start modelling a simplified pH equilibrium, at the interspatial level (right part of section 
(c) in Figure 2), without considering H2O, and 0H~ \ therefore we tuned the reaction rates in order 
to reach the correct percentages of NH^, over total NH^ /NH^ in the different compartments. Like these 
all the rates and initial terms used in this work are obtained by manual adjustments made looking at the 
simulations results and trying to keep simulations times acceptable - we plan to refine these rates and 
numbers in future works to reflect more deeply the available biological data. Following [18], we consider 
an extracellular pH of 4.5 [19]. In such conditions, the percentage of molecules of NH^ over the sum 
NH3 +NH^ should be around 0.002. The reaction we considered is the following: 

with ki = 0.018 * 10"^ and ^2 = 0.562 * 10"^. One can translate this reaction with the SCLS rules. 

NH^hNH^ (Rl) 

NH^^NHs (R2) 

In Figure 3 we show the results of this first simulation given the initial term T = NHi, x 138238 \NH^ x 
138238. 

This equilibrium is different at the intracellular level (pH around 7 and 8) [15], so we use two new 
rules to model the transformations of NH^ and NH^ inside the cell, namely: 

NH3 Jnside NH^ Jnside (R3) 



where = 0.562* 10-^. 



A^i// Jnside A^i/s Jnside 



(R4) 





3.2 LjAMT2;2 Uptake 

We can now present the SCLS model of the uptake of the LjAMT2;2 transporter (left part of section (c) in 
Figure 2). We add a looping sequence modelling an arbusculated plant cell. Since we are only interested 
in the work done by the LjAMT2;2 transporter, we consider a membrane containing this single element. 
The work of the transporter is modelled by the rule: 

NH^ I (L jAMTlf J (X) 3V //+ I (L jAMTlf J {X \ NH3 Jnside) (R5) 
whereyfc^ = 0.1*10-^ 

We can investigate the uptake rate of the transporter at different initial concentrations of NH^ and 
NH^. Figure 4 and Figure 5 show the results for the initial terms 

Ti = NH3 X 776 1 NH^ X 276400 1 (L jAMTlf J e and T2 = NH3 x 276400 1 NH^ x 776 1 (L jAMTlf J e, 
respectively (the left one represents the whole simulations, while on the right there is a magnification of 
their initial segment). 

We can also investigate the uptake rate of the transporter at different extracellular pH. Namely, we 
consider an extracellular pH equal to the intracellular one (pH around 7 and 8), obtained by imposing 
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Figure 6: At extracellular pH=7. 



Rl and R2 equal to /?3 and 7^4, respectively, i.e. ^2 = ^2- Figure 6 shows the results for the initial terms 
Ti =NH3 X 138238 I A^i^+ x 138238 | {LjAMT2f\e. 

Since now we modeled the transporter supposing that no active form of energy is required to do 
the actual work - which means that the NH^ gradient between the cell and the extracellular ambient is 
sufficient to determine a net uptake. The predicted tridimensional structure of LjAMT2;2 suggests that 
it does not use ATP ^ as an energy source [18], nevertheless trying to model an "energy consumption" 
scenario is interesting to make some comparisons. Since this is only a proof of concept there is no need 
to specify here in which form this energy is going to be provided, furthermore as long as we are only 
interested in comparing the initial rates of uptake we can avoid defining rules that regenerate energy in 
the cell. Therefore, rule R5 modelling the transporter role can be modified as follows: 

NH^ I {LjAMT2f\ {ENERGY \X) ^ \ {LjAMTl)^ \ (X |7V//3_inside) (R5') 

which consumes an element of energy within the cell. We also make this reaction slower, since it is 
now catalysed by the concentration of the ENERGY element, actually, we set kff = 0.1 * 10~^^. Given 
the initial term T = NH3 x 138238 |A^7/+ x 138238 1 {LjAMT2f \ {ENERGY x 100000) we obtain the 
simulation result in Figure 7. Note that the uptake work of the transporter terminates when the ENERGY 
inside the cell is completely exhausted. 

3.3 NH^ Diffusing from the Fungus 

We now model the diffusion of NH^ from the fungus to the extracellular level (sections (al), and (a2) 
of Figure 2). In section (al) of the figure, the passage of NH^ to the interfacial periarbuscular space 
happens by diffusion. We can model this phenomenon by adding a new compartment, representing the 
fungus, from which NH^ flows towards the fungus-plant interface. This could be modelled through the 
rule: 

{EungMembrf \ {NH^ | X) NH^ \ {EungMembrf \ X (R6) 

By varying the value of the rate kf one might model different extemalization speeds and thus test 
different hypotheses about the underlying mechanism. In Figure 8 we give the simulation result, with 

^ATP is the "molecular unit of currency" of intracellular energy transfer [26] and is used by many transporters that work 
against chemical gradients. 
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Figure 7: LjAMT2;2 with active energy. 
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Figure 8: Diffusing NH^ from the fungus, kf = I. 
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three different magnification levels, going from the whole simulation on the left to the very first seconds 
on the right, obtained from the initial term Tf (FungMembr)^ \ {NH^ x 2764677) | {LjAMTl)^ \ e 
with kf = I. In the initial part, one can see how fast, in this case, NH^ diffuses into the periarbuscular 
space (in the figure, NH4 represents the quantity of NH^ in that part of the system). In Figure 9 we 
give the simulation result obtained from the same initial term Tf with a slower diffusion rate, namely 
kf = 0.01 * 10~^ (note that the magnification levels in the three panels are different with respect to those 
in Figure 8). 

Finally, we would like to remark, without going into the simulation details, how we can model in 
a rather natural way the portion (a2) of Figure 2 in SCLS. Namely, we need some rules to produce 
vesicle containing NH^ molecules within the fungal cell. Once the vesicle is formed, another rule drives 
its exocytosis towards the interfacial space, and thus the diffusion of the previously encapsulated NH^ 
molecules. The needed rules are given in the following: 



{FungMembrf \ {NH^\X) h {FungMembrf \ (X \ {Vesiclef \NH^) 
{FungMembrf\ (NH^ \ Y \ (Vesicle)^ \X) h {FungMembrf \ {Y \ (Vesicle)^] (NH^ \X)) 



(R7) 

(R8) 
(R9) 



{FungMembry \ {Y \ {Vesicley \X) ^ X \ {FungMembry \ Y 

Where rule R7 models the creation of a vesicle, rule R8 model the encapsulation of an NH^ molecule 
within the vesicle and rule R9 models the exocytosis of the vesicle content. 
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4 Discussion 



We dissected the route for the passage of NH^ I NH^ from the fungus to the plant in known and hy- 
pothetical mechanisms which were transformed in rules. Further, also the properties of the different 
compartments and their influence on the transported molecules were included, thus giving a first model 
for the simulation of the nutrients transfer. With the model so far we can simulate the behaviour of the 
system when varying parameters as the different compartments pH, the initial substrate concentrations, 
the transport/diffusion speeds and the energy supply. 

We can start comparing the two simulations with the plant cell with the Lj AMT2;2 transporter placed 
in different extracellular situations: low NH^ concentration (Figure 5) or high NH^ concentration (Fig- 
ure 4). As a natural consequence of the greater concentration, the ammonium uptake is faster when 
the simulation starts with more NH^, as long as the LjAMT2;2 can readily import it. The real situa- 
tion should be similar to this simulation, assuming that the level of extracellular NH^ I NH^ is stable, 
meaning an active symbiosis. 

The simulation which represents an extracellular pH around 7 (Figure 6) shows a decreased inter- 
nalisation speed with respect to the simulation in Figure 4, as could be inferred from the concentrations 
of NH4 Jnside and NH3_inside in the plots on the right (focusing on the initial activity): this supports 
experimental data about the pH-dependent activity of the transporter and suggests that the extracellular 
pH is fundamental to achieve a sufficient ammonium uptake for the plants. It's noticeable how the initial 
uptake rate in this case is higher, despite the neutral pH, than the rate obtained considering an "energy 
quantum" used by the transporter (which has the same starting term), as could be seen in the right panels 
of Figure 7 and Figure 6. These results could enforce the biological hypothesis that, instead of ATP, a 
NH^ concentration gradient (possibly created by the fungus) is used as energy source by the LjAMT2;2 
protein. 

The simulations which also consider the fungal counterpart are interesting because they provide an 
initial investigation of this rather poorly characterized side of the symbiosis and confirm that plants can 
efficiently gain ammonium if NH^ is released from the fungi. This evidence supports the last biological 
hypothesis about how fungi supply nitrogen to plants [16, 10], and could lead to further models which 
could suggest which is the needed rate for NH^ transport from fungi to the interfacial apoplast; thus 
driving biologists toward one (or some) of the nowadays considered hypotheses (active transport of 
NH^, vesicle formation, etc.). 
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5 Conclusions and Future Work 

This paper reports on the use of SCSL to simulate and understand some biological behaviour which is 
still unclear to biologists, and it also constitues a first attempt to predict biological behaviour and give 
some directions to biologists for future experiments. 

SCLS has been particularly suitable to model the symbiosis (where substances flow through different 
cells) thanks to its feature to model compartments, and their membranes, in a simple and natural way. Our 
simulations have confirmed some of the latest experimental results about the LjAMT2;2 transporter [18] 
and also support some of the hypotheses about the energy source for the transport. These are the first 
steps towards a complete simulation of the symbiosis and open some interesting paths that could be 
followed to better understand the nutrient exchange. 

As demonstrated by heterologous complementation experiments in yeast, mycorrhiza-specific plant 
transporters [20, 18] show different uptake efficiencies under varying pH conditions. Looking on this and 
the results from the simulations with different pH conditions in the periarbuscular space new experiments 
for a determination of the uptake kinetics (TT^-values) under a range of pH seems mandatory. Another 
conclusion from the model is that, for accurate simulation, exact in vivo concentration measurements 
have to be carried out even though at the moment this is a difficult task due to technical limitations. 

As shown by various studies, many transporters on the plant side show a strong transcriptionally 
regulation and the majority of them are thought to be localized at the plant- fungus interface [17]. Con- 
sequently a quantification of these proteins in the membrane will be a prerequisite for an accurate future 
model. 

It is known that some of these transporters (e.g. PT4 phosphate transporters) obtain energy from 
proton gradients established by proton pumps (Figure 2) whereas others as the LjAMT2;2 and aqua- 
porines [14] use unknown energy supplies or are simply facilitators of diffusion events along gradients. 
Thus they conserve the membranes electrochemical potential for the before mentioned proton dependent 
transport processes. To make the story even more complex some of these transporters which are known 
to be regulated in the AM symbiosis (e.g. aquaporines) show overlaps in the selectivity of their sub- 
strates [23, 24, 43, 31]. Consequently integration of interrelated transporters in a future model will be a 
necessary and challenging task. It is quite probable that also transporters for other macronutrients (potas- 
sium and sulfate) which might be localized in the periarbuscular membrane influence the electrochemical 
gradients for the known transport processes. 

With the ongoing sequencing work [27] on the arbuscular model fungus Glomus intraradices trans- 
porters on the fungal side of the symbiotic compartment are likely to be identified and characterized 
soon. Data from such future experiments could be integrated in the model and help to answer the question 
whether transporter mediated diffusion or vesicle based excretion events lead to the release of ammonium 
into the periarbuscular space [10]. 

Further questions about the plant nutrient uptake and competitive fungal reimport process [10, 5] 
might be answered. Based on the transport properties of orthologous transporters from different fungal 
and/or plant species, theories could be developed which explain different mycorrhiza responsiveness of 
host plants; meaning why certain plant- AM fungus combinations have a rather disadvantageous than 
beneficial effect for the plant. 

In future research on AM and membrane transport processes in general values from measurements 
of concentrations or kinetics can and have to be included in the model and will show how the whole 
system is influenced by these values. Vice versa simulations could be the base for new hypotheses and 
experiments. As a future extension of the modelling technique, we plan to follow the direction taken in 
[12, 4] using type systems to guarantee that an SCLS term satisfies certain biological properties (enforced 
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by a set of typing rules), and also investigate how type systems could enrich the study of quantitative 
systems. 
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